(* Content-type: application/vnd.wolfram.mathematica *)

(*** Wolfram Notebook File ***)
(* http://www.wolfram.com/nb *)

(* CreatedBy='Mathematica 8.0' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       157,          7]
NotebookDataLength[     29060,        872]
NotebookOptionsPosition[     26748,        787]
NotebookOutlinePosition[     27087,        802]
CellTagsIndexPosition[     27044,        799]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{

Cell[CellGroupData[{
Cell[BoxData[{
 RowBox[{
  RowBox[{"A", "=", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
       RowBox[{"-", "2.1978"}], ",", "0", ",", "2.1978", ",", 
       RowBox[{"-", "0.6593"}], ",", " ", "0", ",", "0.6593"}], "}"}], ",", 
     RowBox[{"{", "\n", 
      RowBox[{"0", ",", 
       RowBox[{"-", "0.7692"}], ",", 
       RowBox[{"-", "0.7692"}], ",", " ", "0.7692", ",", "0.7692", ",", "0"}],
       "}"}]}], "}"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{"A", "//", "MatrixForm"}]}], "Input",
 CellChangeTimes->{{3.546763699016968*^9, 3.5467637733152*^9}, {
  3.546763823877201*^9, 3.5467638306757717`*^9}}],

Cell[BoxData[
 TagBox[
  RowBox[{"(", "\[NoBreak]", GridBox[{
     {
      RowBox[{"-", "2.1978`"}], "0", "2.1978`", 
      RowBox[{"-", "0.6593`"}], "0", "0.6593`"},
     {"0", 
      RowBox[{"-", "0.7692`"}], 
      RowBox[{"-", "0.7692`"}], "0.7692`", "0.7692`", "0"}
    },
    GridBoxAlignment->{
     "Columns" -> {{Center}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, 
      "RowsIndexed" -> {}},
    GridBoxSpacings->{"Columns" -> {
        Offset[0.27999999999999997`], {
         Offset[0.7]}, 
        Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> {
        Offset[0.2], {
         Offset[0.4]}, 
        Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}],
  Function[BoxForm`e$, 
   MatrixForm[BoxForm`e$]]]], "Output",
 CellChangeTimes->{3.546763831263091*^9, 3.546764591557541*^9, 
  3.546764737986205*^9, 3.5467777171096277`*^9, 3.546778138227702*^9, 
  3.5467785926203003`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"H", "=", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
       RowBox[{"-", "1"}], ",", "0"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"0", ",", 
       RowBox[{"-", "1"}]}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"0", ",", "0"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"0", ",", "0"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"0", ",", "0"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"0", ",", "0"}], "}"}]}], "}"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{"H", "//", "MatrixForm"}]}], "Input",
 CellChangeTimes->{{3.54676378919948*^9, 3.546763839986713*^9}}],

Cell[BoxData[
 TagBox[
  RowBox[{"(", "\[NoBreak]", GridBox[{
     {
      RowBox[{"-", "1"}], "0"},
     {"0", 
      RowBox[{"-", "1"}]},
     {"0", "0"},
     {"0", "0"},
     {"0", "0"},
     {"0", "0"}
    },
    GridBoxAlignment->{
     "Columns" -> {{Center}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, 
      "RowsIndexed" -> {}},
    GridBoxSpacings->{"Columns" -> {
        Offset[0.27999999999999997`], {
         Offset[0.7]}, 
        Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> {
        Offset[0.2], {
         Offset[0.4]}, 
        Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}],
  Function[BoxForm`e$, 
   MatrixForm[BoxForm`e$]]]], "Output",
 CellChangeTimes->{3.5467638404321117`*^9, 3.5467645917450733`*^9, 
  3.546764738098176*^9, 3.546777717261932*^9, 3.546778138323722*^9, 
  3.5467785927162046`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"d", "=", 
   RowBox[{"{", 
    RowBox[{"0", ",", "0", ",", "0.46", ",", "0", ",", "0.46", ",", 
     RowBox[{"-", "0.138"}]}], "}"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{"d", "//", "MatrixForm"}]}], "Input",
 CellChangeTimes->{{3.546763844710994*^9, 3.5467638674510517`*^9}, {
  3.5467795942805853`*^9, 3.5467795993329287`*^9}}],

Cell[BoxData[
 TagBox[
  RowBox[{"(", "\[NoBreak]", 
   TagBox[GridBox[{
      {"0"},
      {"0"},
      {"0.46`"},
      {"0"},
      {"0.46`"},
      {
       RowBox[{"-", "0.138`"}]}
     },
     GridBoxAlignment->{
      "Columns" -> {{Center}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}},
        "RowsIndexed" -> {}},
     GridBoxSpacings->{"Columns" -> {
         Offset[0.27999999999999997`], {
          Offset[0.5599999999999999]}, 
         Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> {
         Offset[0.2], {
          Offset[0.4]}, 
         Offset[0.2]}, "RowsIndexed" -> {}}],
    Column], "\[NoBreak]", ")"}],
  Function[BoxForm`e$, 
   MatrixForm[BoxForm`e$]]]], "Output",
 CellChangeTimes->{3.546763867790099*^9, 3.546764591852703*^9, 
  3.546764738118105*^9, 3.54677771732276*^9, 3.546778138412825*^9, 
  3.5467785927893963`*^9, 3.54677960055794*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"A", ".", "d"}]], "Input",
 CellChangeTimes->{{3.546777727841982*^9, 3.5467777298837023`*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"0.9200046`", ",", "0.`"}], "}"}]], "Output",
 CellChangeTimes->{3.546777730366817*^9, 3.546778138504305*^9, 
  3.5467785928095407`*^9, 3.546779601612137*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"A", ".", "H"}]], "Input",
 CellChangeTimes->{{3.546777731787745*^9, 3.54677773255612*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{"2.1978`", ",", "0.`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0.`", ",", "0.7692`"}], "}"}]}], "}"}]], "Output",
 CellChangeTimes->{3.54677773317726*^9, 3.546778138593473*^9, 
  3.546778592828861*^9, 3.5467796031529303`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"Dnn", "=", "1"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Dss", "=", "1"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Dhat", "=", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{"Dnn", ",", "0"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"0", ",", "Dss"}], "}"}]}], "}"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{"Dhat", "//", "MatrixForm"}]}], "Input",
 CellChangeTimes->{{3.546763896613258*^9, 3.546763921482533*^9}}],

Cell[BoxData[
 TagBox[
  RowBox[{"(", "\[NoBreak]", GridBox[{
     {"1", "0"},
     {"0", "1"}
    },
    GridBoxAlignment->{
     "Columns" -> {{Center}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, 
      "RowsIndexed" -> {}},
    GridBoxSpacings->{"Columns" -> {
        Offset[0.27999999999999997`], {
         Offset[0.7]}, 
        Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> {
        Offset[0.2], {
         Offset[0.4]}, 
        Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}],
  Function[BoxForm`e$, 
   MatrixForm[BoxForm`e$]]]], "Output",
 CellChangeTimes->{3.546763922261869*^9, 3.5467645918779182`*^9, 
  3.546764738137095*^9, 3.5467777173380537`*^9, 3.5467781386762943`*^9, 
  3.5467785928478518`*^9, 3.546779605866314*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"etilde", "=", 
  RowBox[{
   RowBox[{"(", 
    RowBox[{
     RowBox[{"en", "^", "2"}], "+", 
     RowBox[{
      FractionBox["Dss", "Dnn"], 
      RowBox[{"es", "^", "2"}]}]}], ")"}], "^", 
   RowBox[{"(", 
    RowBox[{"1", "/", "2"}], ")"}]}]}]], "Input",
 CellChangeTimes->{{3.546763934980266*^9, 3.546763997540638*^9}, {
  3.546764048101098*^9, 3.546764075985674*^9}}],

Cell[BoxData[
 SqrtBox[
  RowBox[{
   SuperscriptBox["en", "2"], "+", 
   SuperscriptBox["es", "2"]}]]], "Output",
 CellChangeTimes->{3.546764077296997*^9, 3.54676459189721*^9, 
  3.5467647381556683`*^9, 3.54677771735546*^9, 3.54677813878232*^9, 
  3.546778592866577*^9}]
}, Open  ]],

Cell[BoxData[{
 RowBox[{
  RowBox[{"g0", "=", "0.91"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"ef", "=", "3"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"g", "[", "x_", "]"}], ":=", 
   RowBox[{"g0", " ", 
    RowBox[{"Exp", "[", 
     FractionBox[
      RowBox[{"-", "x"}], "ef"], "]"}]}]}], ";"}]}], "Input",
 CellChangeTimes->{{3.546764081301738*^9, 3.5467641266516447`*^9}, {
  3.5467641815953293`*^9, 3.5467641955187397`*^9}, {3.546778585865238*^9, 
  3.546778589255455*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"F", "[", "x_", "]"}], ":=", 
   RowBox[{"Dnn", " ", 
    FractionBox["x", 
     RowBox[{"g", "[", "x", "]"}]]}]}], ";"}]], "Input",
 CellChangeTimes->{{3.546764204006425*^9, 3.546764220331476*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"\[Gamma]", "=", 
  RowBox[{"F", "[", "etilde", "]"}]}]], "Input",
 CellChangeTimes->{{3.546764250747282*^9, 3.546764262518156*^9}}],

Cell[BoxData[
 RowBox[{"1.0989010989010988`", " ", 
  SuperscriptBox["\[ExponentialE]", 
   RowBox[{
    FractionBox["1", "3"], " ", 
    SqrtBox[
     RowBox[{
      SuperscriptBox["en", "2"], "+", 
      SuperscriptBox["es", "2"]}]]}]], " ", 
  SqrtBox[
   RowBox[{
    SuperscriptBox["en", "2"], "+", 
    SuperscriptBox["es", "2"]}]]}]], "Output",
 CellChangeTimes->{3.546764263309856*^9, 3.546764591995837*^9, 
  3.546764738256406*^9, 3.546777717447999*^9, 3.546778138907988*^9, 
  3.546778592944509*^9}]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  RowBox[{"e", "=", 
   RowBox[{"{", 
    RowBox[{"en", ",", "es"}], "}"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.5467643478074512`*^9, 3.546764359319519*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"eqn", "=", 
  RowBox[{
   RowBox[{
    FractionBox[
     RowBox[{"Dhat", ".", "e"}], "\[Gamma]"], "-", 
    RowBox[{"A", ".", "d"}], "+", 
    RowBox[{"A", ".", "H", ".", "e"}]}], "//", "Simplify"}]}]], "Input",
 CellChangeTimes->{{3.546764232025268*^9, 3.546764275817031*^9}, {
  3.546764338552412*^9, 3.546764390333373*^9}, {3.546778085206892*^9, 
  3.5467781694459667`*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{
    RowBox[{"-", "0.9200046`"}], "+", 
    RowBox[{"en", " ", 
     RowBox[{"(", 
      RowBox[{"2.1978`", "\[VeryThinSpace]", "+", 
       FractionBox[
        RowBox[{"0.9100000000000001`", " ", 
         SuperscriptBox["\[ExponentialE]", 
          RowBox[{
           RowBox[{"-", 
            FractionBox["1", "3"]}], " ", 
           SqrtBox[
            RowBox[{
             SuperscriptBox["en", "2"], "+", 
             SuperscriptBox["es", "2"]}]]}]]}], 
        SqrtBox[
         RowBox[{
          SuperscriptBox["en", "2"], "+", 
          SuperscriptBox["es", "2"]}]]]}], ")"}]}]}], ",", 
   RowBox[{"0.`", "\[VeryThinSpace]", "+", 
    RowBox[{"es", " ", 
     RowBox[{"(", 
      RowBox[{"0.7692`", "\[VeryThinSpace]", "+", 
       FractionBox[
        RowBox[{"0.9100000000000001`", " ", 
         SuperscriptBox["\[ExponentialE]", 
          RowBox[{
           RowBox[{"-", 
            FractionBox["1", "3"]}], " ", 
           SqrtBox[
            RowBox[{
             SuperscriptBox["en", "2"], "+", 
             SuperscriptBox["es", "2"]}]]}]]}], 
        SqrtBox[
         RowBox[{
          SuperscriptBox["en", "2"], "+", 
          SuperscriptBox["es", "2"]}]]]}], ")"}]}]}]}], "}"}]], "Output",
 CellChangeTimes->{{3.546764374135219*^9, 3.5467643909029207`*^9}, 
   3.5467645920286207`*^9, 3.546764738294682*^9, 3.546777717568611*^9, {
   3.5467781086936502`*^9, 3.5467781389431353`*^9}, 3.546778170025362*^9, 
   3.546778593023218*^9, 3.546779617675188*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"eqn", "[", 
  RowBox[{"[", "1", "]"}], "]"}]], "Input",
 CellChangeTimes->{{3.546764704874363*^9, 3.546764707513907*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"-", "0.9200046`"}], "+", 
  RowBox[{"en", " ", 
   RowBox[{"(", 
    RowBox[{"2.1978`", "\[VeryThinSpace]", "+", 
     FractionBox[
      RowBox[{"0.9100000000000001`", " ", 
       SuperscriptBox["\[ExponentialE]", 
        RowBox[{
         RowBox[{"-", 
          FractionBox["1", "3"]}], " ", 
         SqrtBox[
          RowBox[{
           SuperscriptBox["en", "2"], "+", 
           SuperscriptBox["es", "2"]}]]}]]}], 
      SqrtBox[
       RowBox[{
        SuperscriptBox["en", "2"], "+", 
        SuperscriptBox["es", "2"]}]]]}], ")"}]}]}]], "Output",
 CellChangeTimes->{3.5467647079424963`*^9, 3.546764738305954*^9, 
  3.546777717615531*^9, 3.546778138964319*^9, 3.54677821725056*^9, 
  3.546778593088586*^9, 3.546779618902582*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"eqn", "[", 
  RowBox[{"[", "2", "]"}], "]"}]], "Input",
 CellChangeTimes->{{3.546764711546939*^9, 3.54676471429808*^9}}],

Cell[BoxData[
 RowBox[{"0.`", "\[VeryThinSpace]", "+", 
  RowBox[{"es", " ", 
   RowBox[{"(", 
    RowBox[{"0.7692`", "\[VeryThinSpace]", "+", 
     FractionBox[
      RowBox[{"0.9100000000000001`", " ", 
       SuperscriptBox["\[ExponentialE]", 
        RowBox[{
         RowBox[{"-", 
          FractionBox["1", "3"]}], " ", 
         SqrtBox[
          RowBox[{
           SuperscriptBox["en", "2"], "+", 
           SuperscriptBox["es", "2"]}]]}]]}], 
      SqrtBox[
       RowBox[{
        SuperscriptBox["en", "2"], "+", 
        SuperscriptBox["es", "2"]}]]]}], ")"}]}]}]], "Output",
 CellChangeTimes->{{3.5467647148413973`*^9, 3.546764738371216*^9}, 
   3.546777717636998*^9, 3.5467781390233383`*^9, 3.546778220623304*^9, 
   3.546778593108922*^9, 3.546779620422365*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"NSolve", "[", 
  RowBox[{
   RowBox[{
    RowBox[{
     RowBox[{"eqn", "[", 
      RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", "0"}], " ", "&&", 
    RowBox[{
     RowBox[{"eqn", "[", 
      RowBox[{"[", "2", "]"}], "]"}], "\[Equal]", "0"}]}], ",", 
   RowBox[{"{", 
    RowBox[{"en", ",", "es"}], "}"}]}], "]"}]], "Input",
 CellChangeTimes->{{3.54676427836246*^9, 3.5467643025799627`*^9}, {
   3.5467643331900997`*^9, 3.546764334837865*^9}, {3.546764400535858*^9, 
   3.546764444778757*^9}, {3.546764516829564*^9, 3.546764541603517*^9}, 
   3.546764602597178*^9, {3.546764632980279*^9, 3.546764648372901*^9}}],

Cell[BoxData[
 RowBox[{"NSolve", "[", 
  RowBox[{
   RowBox[{
    RowBox[{
     RowBox[{
      RowBox[{"-", "0.9200046`"}], "+", 
      RowBox[{"en", " ", 
       RowBox[{"(", 
        RowBox[{"2.1978`", "\[VeryThinSpace]", "+", 
         FractionBox[
          RowBox[{"0.9100000000000001`", " ", 
           SuperscriptBox["\[ExponentialE]", 
            RowBox[{
             RowBox[{"-", 
              FractionBox["1", "3"]}], " ", 
             SqrtBox[
              RowBox[{
               SuperscriptBox["en", "2"], "+", 
               SuperscriptBox["es", "2"]}]]}]]}], 
          SqrtBox[
           RowBox[{
            SuperscriptBox["en", "2"], "+", 
            SuperscriptBox["es", "2"]}]]]}], ")"}]}]}], "\[Equal]", "0"}], "&&", 
    RowBox[{
     RowBox[{"0.`", "\[VeryThinSpace]", "+", 
      RowBox[{"es", " ", 
       RowBox[{"(", 
        RowBox[{"0.7692`", "\[VeryThinSpace]", "+", 
         FractionBox[
          RowBox[{"0.9100000000000001`", " ", 
           SuperscriptBox["\[ExponentialE]", 
            RowBox[{
             RowBox[{"-", 
              FractionBox["1", "3"]}], " ", 
             SqrtBox[
              RowBox[{
               SuperscriptBox["en", "2"], "+", 
               SuperscriptBox["es", "2"]}]]}]]}], 
          SqrtBox[
           RowBox[{
            SuperscriptBox["en", "2"], "+", 
            SuperscriptBox["es", "2"]}]]]}], ")"}]}]}], "\[Equal]", "0"}]}], 
   ",", 
   RowBox[{"{", 
    RowBox[{"en", ",", "es"}], "}"}]}], "]"}]], "Output",
 CellChangeTimes->{{3.546764280443454*^9, 3.546764292997815*^9}, 
   3.546764335272443*^9, 3.546764434313555*^9, 3.546764483643256*^9, {
   3.546764530599207*^9, 3.546764542878292*^9}, {3.546764593019991*^9, 
   3.546764604616187*^9}, {3.546764641220497*^9, 3.546764649670794*^9}, 
   3.546764739366118*^9, 3.54677771870252*^9, 3.54677813949014*^9, 
   3.546778222689437*^9, 3.546778593695052*^9, 3.546779629222019*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"NSolve", "[", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{
     RowBox[{
      RowBox[{
       RowBox[{"x", "^", "5"}], "-", 
       RowBox[{"2", "x"}], "+", "3"}], "\[Equal]", "0"}], ",", 
     RowBox[{
      RowBox[{"x", "+", "y"}], "\[Equal]", "1"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"x", ",", "y"}], "}"}]}], "]"}]], "Input",
 CellChangeTimes->{{3.546764505921891*^9, 3.5467645101664143`*^9}, {
  3.546764555893607*^9, 3.5467645749141293`*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{
     RowBox[{"x", "\[Rule]", 
      RowBox[{"-", "1.4236058485523289`"}]}], ",", 
     RowBox[{"y", "\[Rule]", "2.423605848552329`"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"x", "\[Rule]", 
      RowBox[{
       RowBox[{"-", "0.2467292569105622`"}], "-", 
       RowBox[{"1.3208163474502472`", " ", "\[ImaginaryI]"}]}]}], ",", 
     RowBox[{"y", "\[Rule]", 
      RowBox[{"1.2467292569105621`", "\[VeryThinSpace]", "+", 
       RowBox[{"1.3208163474502472`", " ", "\[ImaginaryI]"}]}]}]}], "}"}], 
   ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"x", "\[Rule]", 
      RowBox[{
       RowBox[{"-", "0.2467292569105622`"}], "+", 
       RowBox[{"1.3208163474502472`", " ", "\[ImaginaryI]"}]}]}], ",", 
     RowBox[{"y", "\[Rule]", 
      RowBox[{"1.2467292569105621`", "\[VeryThinSpace]", "-", 
       RowBox[{"1.3208163474502472`", " ", "\[ImaginaryI]"}]}]}]}], "}"}], 
   ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"x", "\[Rule]", 
      RowBox[{"0.9585321811867296`", "\[VeryThinSpace]", "-", 
       RowBox[{"0.4984277790318463`", " ", "\[ImaginaryI]"}]}]}], ",", 
     RowBox[{"y", "\[Rule]", 
      RowBox[{"0.04146781881327044`", "\[VeryThinSpace]", "+", 
       RowBox[{"0.4984277790318463`", " ", "\[ImaginaryI]"}]}]}]}], "}"}], 
   ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"x", "\[Rule]", 
      RowBox[{"0.9585321811867296`", "\[VeryThinSpace]", "+", 
       RowBox[{"0.4984277790318463`", " ", "\[ImaginaryI]"}]}]}], ",", 
     RowBox[{"y", "\[Rule]", 
      RowBox[{"0.04146781881327044`", "\[VeryThinSpace]", "-", 
       RowBox[{"0.4984277790318463`", " ", "\[ImaginaryI]"}]}]}]}], "}"}]}], 
  "}"}]], "Output",
 CellChangeTimes->{
  3.5467645105186377`*^9, {3.546764557125167*^9, 3.546764593113366*^9}, 
   3.5467647394775667`*^9, 3.5467777188658752`*^9, 3.546778139590114*^9, 
   3.546778225039752*^9, 3.546778593802299*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"eqn", "[", 
    RowBox[{"[", "2", "]"}], "]"}], "/.", 
   RowBox[{"{", 
    RowBox[{"es", "\[Rule]", "0"}], "}"}]}], "//", "Simplify"}]], "Input",
 CellChangeTimes->{{3.546764755357942*^9, 3.546764774043178*^9}}],

Cell[BoxData["0.`"], "Output",
 CellChangeTimes->{{3.546764767486311*^9, 3.5467647743751183`*^9}, 
   3.546777719183567*^9, 3.546778139672114*^9, 3.5467782262591057`*^9, 
   3.546778594021783*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"eqn", "[", 
    RowBox[{"[", "1", "]"}], "]"}], "/.", 
   RowBox[{"{", 
    RowBox[{"es", "\[Rule]", "0"}], "}"}]}], "//", "Simplify"}]], "Input",
 CellChangeTimes->{{3.546764782218669*^9, 3.546764782357521*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"-", "0.9200046`"}], "+", 
  RowBox[{"en", " ", 
   RowBox[{"(", 
    RowBox[{"2.1978`", "\[VeryThinSpace]", "+", 
     FractionBox[
      RowBox[{"0.9100000000000001`", " ", 
       SuperscriptBox["\[ExponentialE]", 
        RowBox[{"-", 
         FractionBox[
          SqrtBox[
           SuperscriptBox["en", "2"]], "3"]}]]}], 
      SqrtBox[
       SuperscriptBox["en", "2"]]]}], ")"}]}]}]], "Output",
 CellChangeTimes->{3.546764782862822*^9, 3.546777719350192*^9, 
  3.546778139754444*^9, 3.546778227380373*^9, 3.546778594107348*^9, 
  3.546779634709736*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Plot", "[", 
  RowBox[{
   RowBox[{
    RowBox[{"-", "0.92"}], "+", " ", 
    RowBox[{"2.1978", " ", "en"}], "+", " ", 
    RowBox[{"0.91", 
     RowBox[{"Exp", "[", 
      RowBox[{"-", 
       FractionBox["en", "3"]}], "]"}]}]}], ",", 
   RowBox[{"{", 
    RowBox[{"en", ",", "0.005", ",", "0.006"}], "}"}]}], "]"}]], "Input",
 CellChangeTimes->{{3.5467648925203667`*^9, 3.546764994215266*^9}, {
  3.546765211234015*^9, 3.5467652178469677`*^9}, {3.546778238045603*^9, 
  3.546778273128171*^9}, {3.5467784173733187`*^9, 3.546778417578998*^9}, {
  3.546778483242898*^9, 3.546778500142118*^9}, {3.5467786002505074`*^9, 
  3.5467786011452103`*^9}, {3.5467796373393707`*^9, 3.546779652184901*^9}, {
  3.546779921196451*^9, 3.546779934985856*^9}}],

Cell[BoxData[
 GraphicsBox[{{}, {}, 
   {Hue[0.67, 0.6, 0.6], LineBox[CompressedData["
1:eJwV0H0w1HkcB/A9bZ668ZiIdJ4GF2VxymzWx97uiTzsUvu7OvKUh2KPLlZa
6+HcuZYja6O4Uk513XVCzhxa9VsnWrQJV4Zj1BzOQ7pV1FyY7vP74zPfec17
5jvvz8c2Li08QYtGo4XgUO8tS7PA5+IMoDnYqPqATXqrXWs8s9FC/m/zTDZ5
P4+7VEg5Ju+y/m42OTaVXr1DguaOZux3ZZPajUNz6Tloeo5tnwWbjOHIirXz
0fHfinuX/EhT4aaHLoVomdxFddWPPH33A15mGbq2vqhbx4+capZLe66jDzFa
tVJYZNxePe9Ld9DKts95r5mk+jxvn94Amp9Cp0m9yREus9h9Gi0VPvCY9SKX
C/U5k/+hNZxypsiTbNqa5xZgJAIa4+mx+2oGeSi/KzvAEc2vkscnuJKbxYus
R0y07K1oXupEpvUoPA3D0E2i6AyOHTlxPe3KkwS0cnWB32hJeh96newpRtPe
JAbqGZPuR6ySg2XofKaJ02fv7vmuOE/OX6PyrDH9i3O+85NjZ7gKKpeK9Co2
QUFNpMxrEP0setTczQJ2vewyqJ9Ga2iOKRobyCofUE6sUn0TKvw9HOHpA+Vs
hVEm9mHVnbJ0gWE+wdVxQlu4Zg95MaCjXvHeiIV2fm+0ccgD2i9qu9eGoasc
hHFsLzg4sVI5n4R+ZrnCmt4DG5r7TtyRoGfNfLmRTLhyOO1X23K0bO+xmR4f
GD3QtpV7A82/3KedCOD0fGZ9QYE2Ornk/MIPTDvbzXhD6A/XakrhU3iRHlka
OEvl93yUJRwYbZ8S966hJWx3aRcX4sONL9FNT2Hf/KCyzf6wIdfes9sJzfir
zYW3D8okG3dY+6BrHUzPfxMA7ervDzLC0ZKIKhgKhLe33y6rEtHKH1Ny9YNg
eFand1sOesoydS0yGLqa103elKNLVhWG1SFQy+vqSPwJrate1jwKBSvY+YO8
Az1Ovy015MPNTskiZxCdqS4+Ws+HkfMnMpr+Qf+poUURYbDk+0lEyzqaH/Ul
rIfBR+KJVpZJFt4rUrxQFQ4LOi4BBc7o1FjHVvYB6HeOCg1moU2ufRUyeQAu
bIttuRWGHitu6jt9ELSO6ob2JaFVxuVJtgL42qfyiFCC3s0zGLgrgJMr+nkq
yi1kapFSAPEN+bUOOWht9wFOlwD8bYR/j1NmbJG1qQSgT+ccD8lDn503rhsW
gLxfk7GrAC1UmInmBHDti+ASzRl0z7q1lRkBlaZ/3AqWonWLcp+YE/Ddwz0D
P1PebjZZZknAMV97k9gi9JBbLd2GgJ127y4MFqN90+xe7iDg97kbV5tL0X5b
nDr9CLhRZ91tcBY91XAmm0NAVcS5mWTKUQGzXv4EiNW5H9uXUf/l/HIziAC4
LWg8J6P2f+VSQRDASO5//C/l7fKS0MME2NmzXwWVo/d7LOpGEkCvdPWiy9F2
ooacWAJWQuqIGMqPzA294wmY0bbI6qAcrkh7lUjACFlabXEOXR/9uP44Aaqs
DYoMyuN09yQhAf8DP8kE3A==
     "]]}},
  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
  Axes->True,
  AxesOrigin->{0.005, 0},
  PlotRange->{{0.005, 0.006}, {-0.000526403440972989, 0.0013686187485981716`}},
  PlotRangeClipping->True,
  PlotRangePadding->{
    Scaled[0.02], 
    Scaled[0.02]}]], "Output",
 CellChangeTimes->{{3.546779923004689*^9, 3.54677993537776*^9}}],

Cell[BoxData[
 GraphicsBox[{{}, {}, 
   {Hue[0.67, 0.6, 0.6], LineBox[CompressedData["
1:eJwVkn081AcAxn+ZUF2xuuVcDPvY1JLuUxHGHgofjVohbzEW2yELG5K8HJ2X
inIScl3pKMdu1528/bRzp+OnSUmWENJKCWUSqZTZH8/n+e/7ff54jPdHuv+k
RhDEzoX838JRMuZWpvG39vrd14cSuIpwtc+UVXQCxJNYj3B/rmKGkSA9RdcB
odc5NGXHVdiWidcm042g7CQnNRZxFQ+ft4TE0llQyj3HzTKPKmTZd9Wi6fZQ
hlc5xp9OUzAHezoP0neDc0A/jSbmKG5r+JO/0INAWLoGeTUcVtAf7fTNaYoC
wTUIHIvYr4jexH+8eDgeRPIZqrrbB5/HnHjREckB5/kybotnLN5F/zgvVnGh
1GbtsRxMg2TdCgb/VhaUe4RHPs84Dm/1Lmc985wFno2lWRgP/Czd1XR5LpTp
HnFaiWfwZoxFM2nMg1LBE464nMXtIK/liaP5MFrZYDTLFoC36Rmvf7IArMJ+
E25JCaIEqujFE0WQZg7yQ48K0aYn9o/T5YMjdx1gnyqDztR82wRTAOXp9H0h
8ZfRORDoKNtxAVE6DyyTj4twX5qebVBXAoKIDJ3Oq4StMa/+ae1FBAU0i3+O
FyOg+ItRj3YhdGyHk2czJJhuCgtw6ihF6H6XYAdrKVxfWXb3PSuDj0WVbtF2
GcwnKhI9318Cwyb+nv6OKmQU3w1xn7sMovS4lb7dVZzbler+eqUIWg7jzRKn
aqyVDacydStgRVOaXrOtwbahiUGT1ZUgjKxmpsxrocfgOHR9/TsI34dHmizq
wP7YMrx4gxgi07hdmuvr8dsL3ql+xz8QNWtt+HgLCVrFJ+xCbwl6lIK/x6wb
wFy1ntPtfQUiabvjgQ3XMLUtoqPfUoqSJ3kuaRZ/Ir9LkJJbL4VOyLy1tpkc
id8JnWu3yhB6pTG21aARqpYwHUmTDPae1R0MEwWcGLT37G1VoKsyXSJpSnwZ
GT0ydasKr/0S/Ax5StgLZDMOu6/CvzLCesnyJghlzWvU7lxFbm7dwayMJgS3
FW/P96uG+szt4r1LryOwwTV75EE1RF9l3NmXfh0NEsX9FaE1qO4KvCBdqoL/
MQ2b8skaiLxLCH6OCtbfrJUv+rUWzT29pXLNZsRM0rabf6gFK/hwxcbsZowm
tz/XT6kDy9DqRK9mC8pf/VCe9LEOPiGizYzUFqjCq9nSrHqwOKFq5eoUmA1P
XC4tISFmZoxFpFD4kKTOHF5GgisoMuWlUhhyMBkzWUFijjMaXHOUwuW/gnPK
VpIgwtwGPmRS2Nz7T2fpGhKseq2uk7kU3GaH/IQbSHDWBbdJL1JIshyIuLCH
RNCy8yOvVRSC3s3ZPfQgIXI6Z6ZHUXBs1Nc29CJR1NcbZXeDwlJnf9l5PxLK
VavepbdTKNj7YFoQTGJcnqSte4+CJKY3+VwsiS03NjtbP6XAs3r7ff8hEv8W
mJ4MGKEQM8cw1k8g4S/3604dpWDD9VHxk0nk66Sw215SaD19X5OfScL+piJ7
3xsKld5vevqOLexzl3SnvKVwco1uJTN7wb+x36jsPQXPUi+34twF3vh0zfg8
ha3sOIO+PBKONx8Rn6q1grm+4KXeGRJl8zQ3C/VWfHhZo/AtJOF27FChr0Yr
hqru5Z49u/DLUOPHSVqt+A/ypFaO
     "]]}},
  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
  Axes->True,
  AxesOrigin->{0, 0},
  PlotRange->{{0, 0.1}, {-0.009999996133741473, 0.17994664755207201`}},
  PlotRangeClipping->True,
  PlotRangePadding->{
    Scaled[0.02], 
    Scaled[0.02]}]], "Output",
 GeneratedCell->False,
 CellAutoOverwrite->False,
 CellChangeTimes->{{3.546764971052758*^9, 3.546764995169757*^9}, {
   3.546765212138378*^9, 3.546765218249819*^9}, 3.546777719489706*^9, 
   3.546778139839315*^9, 3.546778273926159*^9, 3.546778431827118*^9, {
   3.546778483858238*^9, 3.546778501102748*^9}, {3.546778594173896*^9, 
   3.546778601633698*^9}, {3.5467796383591967`*^9, 3.546779652565132*^9}, 
   3.546779917200158*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"NSolve", "[", 
  RowBox[{
   RowBox[{
    RowBox[{
     RowBox[{"-", "1"}], "+", " ", 
     RowBox[{"2.1978", " ", "en"}], "+", " ", 
     RowBox[{"0.91", 
      RowBox[{"Exp", "[", 
       RowBox[{"-", 
        FractionBox["en", "3"]}], "]"}]}]}], "\[Equal]", "0"}], ",", 
   RowBox[{"{", "en", "}"}]}], "]"}]], "Input",
 CellChangeTimes->{{3.5467651925971527`*^9, 3.5467651930915327`*^9}, {
  3.5467782907637243`*^9, 3.546778309145545*^9}, {3.546778506524487*^9, 
  3.546778506697115*^9}, {3.546778607464615*^9, 3.546778608135087*^9}}],

Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"NSolve", "::", "ifun"}], "MessageName"], 
  RowBox[{
  ":", " "}], "\<\"Inverse functions are being used by \[NoBreak]\\!\\(NSolve\
\\)\[NoBreak], so some solutions may not be found; use Reduce for complete \
solution information. \\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", \
ButtonStyle->\\\"Link\\\", ButtonFrame->None, \
ButtonData:>\\\"paclet:ref/NSolve\\\", ButtonNote -> \\\"NSolve::ifun\\\"]\\)\
\"\>"}]], "Message", "MSG",
 CellChangeTimes->{{3.546778304323251*^9, 3.546778309616867*^9}, 
   3.5467785072060003`*^9, {3.546778594346368*^9, 3.546778608739193*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{"en", "\[Rule]", 
     RowBox[{"-", "9.556508640888438`"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"en", "\[Rule]", "0.04744701376177118`"}], "}"}]}], 
  "}"}]], "Output",
 CellChangeTimes->{{3.54677830437129*^9, 3.5467783096241198`*^9}, 
   3.546778507211557*^9, {3.5467785943522778`*^9, 3.5467786087423077`*^9}}]
}, Open  ]]
},
WindowSize->{740, 623},
WindowMargins->{{Automatic, -1207}, {Automatic, 2}},
FrontEndVersion->"8.0 for Linux x86 (64-bit) (February 23, 2011)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[CellGroupData[{
Cell[579, 22, 634, 16, 69, "Input"],
Cell[1216, 40, 923, 24, 57, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[2176, 69, 647, 20, 50, "Input"],
Cell[2826, 91, 861, 26, 121, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[3724, 122, 370, 8, 50, "Input"],
Cell[4097, 132, 895, 27, 121, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[5029, 164, 117, 2, 30, "Input"],
Cell[5149, 168, 197, 4, 30, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[5383, 177, 114, 2, 30, "Input"],
Cell[5500, 181, 299, 8, 30, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[5836, 194, 503, 15, 88, "Input"],
Cell[6342, 211, 773, 20, 57, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[7152, 236, 395, 12, 54, "Input"],
Cell[7550, 250, 271, 7, 41, "Output"]
}, Open  ]],
Cell[7836, 260, 515, 14, 91, "Input"],
Cell[8354, 276, 244, 7, 55, "Input"],
Cell[CellGroupData[{
Cell[8623, 287, 155, 3, 30, "Input"],
Cell[8781, 292, 509, 15, 43, "Output"]
}, Open  ]],
Cell[9305, 310, 186, 5, 30, "Input"],
Cell[CellGroupData[{
Cell[9516, 319, 401, 10, 58, "Input"],
Cell[9920, 331, 1540, 43, 142, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[11497, 379, 145, 3, 30, "Input"],
Cell[11645, 384, 776, 22, 73, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[12458, 411, 144, 3, 30, "Input"],
Cell[12605, 416, 778, 21, 73, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[13420, 442, 634, 15, 30, "Input"],
Cell[14057, 459, 1924, 51, 142, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[16018, 515, 480, 14, 30, "Input"],
Cell[16501, 531, 1947, 49, 107, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[18485, 585, 259, 7, 30, "Input"],
Cell[18747, 594, 196, 3, 30, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[18980, 602, 259, 7, 30, "Input"],
Cell[19242, 611, 598, 17, 85, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[19877, 633, 766, 17, 52, "Input"],
Cell[20646, 652, 1995, 38, 205, "Output"],
Cell[22644, 692, 2476, 47, 223, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[25157, 744, 561, 14, 52, "Input"],
Cell[25721, 760, 620, 12, 61, "Message"],
Cell[26344, 774, 388, 10, 30, "Output"]
}, Open  ]]
}
]
*)

(* End of internal cache information *)
